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Several recent studies of energy transfer in photosynthetic light harvesting complexes have revealed 
a subtle interplay between coherent and decoherent dynamic contributions to the overall transfer 
efficiency in these open quantum systems. In this work we systematically investigate the impact of 
temporal and spatial correlations in environmental fluctuations on excitation transport in the Fenna- 
Matthews-Olson photosynthetic complex. We demonstrate that the exact nature of the correlations 
can have a large impact on the efficiency of light harvesting. In particular, we find that (i) spatial 
correlations can enhance coherences in the site basis while at the same time slowing transport, and 
(ii) the overall efficiency of transport is optimized at a finite temporal correlation that produces 
maximum overlap between the environmental power spectrum and the excitonic energy differences, 
which in turn results in enhanced driving of transitions between excitonic states. 



I. INTRODUCTION 



Recent experimental revelations of long-lived elec- 
tronic coherence in photosynthetic light harvesting sys- 
tems [IH1] and conjugated polymers [5] have prompted 
a renewed examination of energy transport in densely 
packed molecular aggregates. In particular, the effect of 
quantum coherent dynamics, and its interplay with envi- 
ronmental decoherence and dissipation have been closely 
scrutinized lately (e.g. [6lll2j). Combined with earlier 
detailed modeling of excitation transport in photosyn- 
thetic light harvesting (e.g. [P3HT5] ). a preliminary un- 
derstanding of the complex dynamics of excitation trans- 
port is being molded. However, the picture is far from 
complete. Most photosynthetic light harvesting com- 
plexes (LHCs) are surrounded by protein structures that 
serve multiple functions, including maintenance of struc- 
tural stability |16j and creation of energy landscapes that 
facilitate energy transfer to reaction centers in the core 
chromophoric networks of LHCs 17]. These same pro- 
tein structures also provide a dynamic environment that 
interacts with the chromophore molecules that carry the 
excitation energy. The dynamics of this environment are 
complex and not very well characterized. Environmental 
fluctuations are generally correlated in time [HI HH US] 
and are also believed to be spatially correlated [J [SJ. 
The effects of these correlations on photosynthetic energy 
transport are generally not well understood (although 
recent experimental results suggest that spatial correla- 
tions may be directly responsible for long-lived electronic 
coherence [SJ[5]). Further, new dynamical models capa- 
ble of simulating some of these correlations (e.g. [2TJ1 |2"T] ) 
suggest that their effects are quite significant. Thus it 
is important to determine the range of possible conse- 
quences of temporal and spatial correlations of environ- 
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mental fluctuations on excitation energy transfer both 
within a single LHC and between networked photosyn- 
thetic units. 

In this paper wc undertake a numerical study of 
these effects, analyzing how the spatiotemporal correla- 
tions of environmental fluctuations affect energy trans- 
port in a single LHC. We use a specific system, the 
Fenna-Matthews-Olson (FMO) bacteriochlorophyll com- 
plex [33] as a prototypical model for our study be- 
cause its structural and energetic properties are particu- 
larly well characterized [2"3l[2"lj . We examine the effects of 
environmental fluctuations in a systematic fashion, em- 
ploying model forms of correlations in time and space in 
order to assess the generic effects of these correlations on 
energy transfer efficiency. The results obtained here with 
FMO can thereby be taken as indicative of the generic 
effects of environmental fluctuations on other light har- 
vesting complexes. 

The FMO complex is a small protein in green sulfur 
bacterium that acts as a highly efficient energy transfer 
wire connecting chlorosomes, i.e., light collecting pigment 
arrays, to photosynthetic reaction centers. Structurally, 
the FMO protein is a trimer whose monomers are be- 
lieved to function independently |25j . Each monomer 
contains seven bacteriochlorophyll- a (BChla) molecules 
embedded within a protein scaffold. Recent studies have 
determined the orientation of the FMO complex within 
the inter-mcmbranc region between the chlorosome an- 
tenna and reaction center [251 HE] ■ They present strong 
evidence that the reaction center is strongly coupled to 
BChla 3 and that the excitation energy enters an FMO 
monomer from the chlorosomes via BChla 1 or BChla 6. 
Below, we shall refer to the individual BChla molecules 
as sites or chromophores. 

Under moderate laser driving (or in vivo) there is 
at most one excitation in a single FMO complex. In 
these conditions, the reversible dynamics of the elec- 
tronic degrees of freedom are described by the so-called 
Frenkel exciton Hamiltonian: H c \ = Ylj=i-Ej\j)(j\ + 
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Ej=iEl>j Jij(\j}(i\ + NXiD- Here |j) represents the 
state where only the jth chromophore is excited and all 
other chromophores are in their electronic ground states. 
Ej is the transition energy of chromophore j, including 
any static shifts due to the interactions with the pro- 
tein environment, and Jij describe the excitonic cou- 
pling between chromophores i and j. We adopt a clas- 
sical stochastic model of the FMO environment and de- 
scribe the interactions between excitations and surround- 
ing protein environments as fluctuations of the chro- 
mophore transition energies: H s = $Zj=i Aj'(*)li)01) 
where Aj (t) are time dependent random variables whose 
properties we shall describe shortly. We ignore fluctua- 
tions in the off-diagonal couplings for simplicity, but our 
treatment below can be generalized to treat such noise as 
well. This model is equivalent to a stochastic Liouville 
treatment of excitation dynamics and will lead to dephas- 
ing in the site basis when the dynamics are averaged over 
the random process [2"7lf2l?] . Since we are treating the 
fluctuating phonon environment as a classical quantity, 
dynamics under this stochastic Liouville treatment are 
only exact for an infinite temperature phonon environ- 
ment, although they are known to become increasingly 
accurate for high temperature baths [29J. In the context 
of this work, the primary utility of this model comes from 
the fact that it allows one to numerically incorporate 
any spatial and temporal correlation of the environment 
into the dynamics of the chromophores, and thus enables 
systematic studies of a range of correlated bath dynam- 
ics. This is difficult or even impossible for alternative 
treatments of the chromophore-environment interactions 
such as Redfield equations, Markovian master equations, 
or cumulant expansion techniques |35| . Since the main 
conclusions we draw from the present study will pertain 
to the relative effects of environmental correlations, this 
stochastic model will be sufficient for our purposes. In 
summary, the dynamics of the FMO complex will be de- 
scribed by the following master equation: 

do -i 7 

dt = T^ H 'd +ll Yj + 7tP[<4ap<^ ]p (1) 

3=1 

where H = H c \ + H s , and the second and third terms 
describe, respectively, radiative excitation decay and ex- 
citation trapping at site 3 due to interaction with the 
reaction center. The Lindblad superoperator V[A]p 
A f) A — ^A^Ap — \pAA A for any operator A, and <rj = 
\0)(j\ is a lowering operator (|0) denotes the electronic 
ground state of all seven chromophores) . 

The stochastic noise processes (i) are Gaussian with 
zero mean, a reasonable form for noise terms resulting 
from coupling to a dynamic protein environment |13j . 
The temporal correlations are taken to be of exponen- 
tial form: (Aj(t)Aj(t + t)) = cxp(-|r|/r c ), where 
the correlation time, r c , is a free parameter which we 
will vary in the simulations below. The magnitude of the 
noise variation, A , characterizes the size of the environ- 
mental fluctuations, which we expect to be a function 



of the environmental dynamics and temperature. Con- 
sistent with the stochastic Liouville treatment, we use 
a high temperature approximation for the variance of 
phononic fluctuations (e.g. see Ch. 8 of Ref. [3D]) to 
arrive at Aq = 2Ej^ksT, where Er is the reorganiza- 
tion energy and T is the temperature. The true form of 
the decay of temporal correlations in protein dynamics is 
complex. It is typically a combination of Gaussian decay 
at initial times and multiple exponential decays at longer 
times. Due to the complexity of simulating such a cor- 
relation and the lack of detailed knowledge of temporal 
correlation in the FMO environment, we have chosen here 
to simulate the temporal correlation decay with a single 
exponential. At the temperatures and timescales relevant 
to light harvesting complex dynamics this approximation 
is reasonable and corresponds to a coarse-graining of en- 
vironmental dynamics [20] . This is also the justification 
provided for the use of the overdamped Brownian oscilla- 
tor model commonly employed to model protein dynam- 
ics in LHCs (e.g. [T31 IHJ E3] ) ■ In fact, the single expo- 
nential decay form we use is the high-temperature limit of 
the symmetrized correlation function for the overdamped 
Brownian oscillator model - where the Matsubara terms 
can be ignored [41 . Another justification for this form 
of temporal correlation decay comes from Doob's theo- 
rem |42j , which states that any stationary, Gaussian and 
Markov random process (all plausible properties for the 
environmental fluctuations in pigment-protein complexes 
at physiological temperatures) must possess an exponen- 
tially decaying temporal correlation. 

Spatial correlations of the noise processes are described 
by a matrix C with elements = -^(Ai(t)Aj(t)). The 
spatial correlation of protein dynamics in photosynthetic 
complexes is not a well studied subject and there are cur- 
rently no precise characterizations of either the nature or 
the extent of spatial correlation between electron-phonon 
couplings or phonon fluctuations. In this study we will 
therefore explore the effects of generic forms for C, within 
the restriction of having positive correlations between 
fluctuations at different chromophores (i.e., positive def- 
inite matrix > 0). Specifically, we will examine the 
following three models: 

1. No spatial correlations. = 

2. Dimerized correlations. Excitons in FMO are mostly 
delocalized on two chromophores [31 j - The chro- 
mophores that are the most strongly coupled are the 
pairs: 1-2, and 5-6, and to a lesser extent, 4-5 and 4- 
7. This dimerization motivates us to use a correlation 
matrix, C D , with the only non-zero off-diagonal entries 
being: Cf 2 = Cg = Cg = Cg = 0.9, Cg = Cg = 
C47 = Cf 4 = 0.4. That is, only the strongly coupled 
chromophores experience correlated fluctuations. 

3. Distance- dependent correlations. The typical dis- 
tances between chromophores in the C. tepidum FMO 
monomer are known from its crystal structure |23l 124] . 
It is common in the literature to assume spatial corre- 
lations that decay exponentially with distance: = 
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FIG. 1: (Color online) Trapping probability at 20 ps as a function of the temporal correlation of the FMO 
environment, for various spatial correlation models. The initial state is an excitation localized on BChla 1, and the 
temperature in both panels is T=77K. The curves are polynomial fits to the data points that indicate the general 

trend, and the error bars show the standard deviation of the average taken over 100 sample evolutions. In both 
panels, the blue (solid) curve corresponds the case of no spatial correlation C", the red (dotted) curve corresponds 
to the case of dimerized spatial correlation C D , and the black (dashed) curves correspond to exponential spatial 
correlation C E , with the labeled correlation radii. Panel (a) compares the trapping probability for the case of no 
spatial correlation to the case of dimerized correlations, and panel (b) compares no correlation with exponential 

correlations. 
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FIG. 2: (Color online) Trapping probability at 20 ps as a function of the temporal correlation of the FMO 
environment, for various spatial correlation models. The initial state is an excitation localized on BChla 1, and the 
temperature in both panels is T=300K. The curves are polynomial fits to the data points that indicate the general 
trend, and the error bars show the standard deviation of the average taken over 100 sample evolutions. The color 
coding for the curves is the same as in Fig. [TJ Panel (a) compares the trapping probability for the case of no spatial 
correlation to the case of dimerized correlations, and panel (b) compares no correlation with exponential correlations. 



e~ di i' Rc , where R c is a variable correlation radius [37]. Taken together, these models 1-3 allow a systematic 
However, microscopic justification for this form is an study of the effects of increasing spatial correlations, 
open question |40| and so we have carried out calcu- 
lations with both exponential and inverse polynomial 
decay of spatial correlations. 
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II. RESULTS 

Eq. ([I]) was numerically integrated using the Adams- 
Bashforth-Moulton multistep method. The temporal 
correlations in the noise processes were generated from 
white noise by passing the latter through a finite im- 
pulse response filter [33] ■ At each time instant the spa- 
tial correlation between fluctuations at different sites was 
realized by the standard technique of using the Cholesky 
decomposition of C, see e.g. [35]. The electronic Hamil- 
tonian, H e \, was formed using site energies and coupling 
strengths for C. tepidum FMO from tables 1 (column 3) 
and 4 (trimer column) of Ref. [35J . Finally, we used val- 
ues of Er = 35cm -1 , ji = 1 ns -1 and 7t = 1 ps -1 , all of 
which are consistent with the most detailed experimental 
and theoretical literature on FMO [25] [32] . 

We focus primarily on the time evolution of the average 
trapping probability of the excitation. This is formally 
defined as: Pt rap (i) = EA[(trap|p(t)|trap)], where pit) 
is the state of the system at time t, and Ea[ - ] denotes 
an ensemble average over instances of the stochastic pro- 
cess. In the simulations presented below, we perform 
this ensemble average over 100 instances. Ptrap provides 
a measure of the efficiency of excitation transport across 
the FMO protein, from the initial absorption at sites 1 
or 6, to the trap, which represents the reaction center. 
Figures [l] - [4j show the average trapping probability at 
20ps as a function of the temporal correlation time of the 
environmental fluctuations for various temperatures and 
initial states. In all the figures, the blue (solid) curve 
shows the trapping probability for the case of no spa- 
tial correlations, the red (dotted) curve shows Ptrap for 
the case of dimerized spatial correlations, and the black 
(dashed) curves shows the same quantity for exponential 
spatial correlations, with the labeled correlation radii. 
We choose to plot the average trapping probability at 
20ps because we wish to compare relative efficiencies at 
times when transient effects are absent and thus are in- 
terested in the asymptotic behavior. However, we note 
that at earlier times the curves look qualitatively similar 
to those in Figs. [T]-[4]and the conclusions given below 
hold for all times in the interval 10 — 25ps. 

Several distinct features are evident from Figs, [I] - 
[4] Firstly, with any spatial correlation, there is an op- 
timal temporal correlation time where the efficiency of 
excitation transport is maximized with respect to envi- 
ronmental fluctuations. For T = UK this optimal time 
is ~ 40fs when the initial state is an excitation on site 1 
(Fig. [I]), and ~ 30fs when the initial state is an excita- 
tion on site 6 (Fig. [3]) . The efficiency dies off very quickly 
for shorter correlation times and more slowly for longer 
correlation times. A second feature evident in these fig- 
ures is the large effect that spatial correlations have on 
the efficiency of excitation transport. Uncorrelated fluc- 
tuations provide the greatest efficiency and generally the 
efficiency decreases with increasing spatial correlation. 
There is very little difference between the uncorrelated 
spatial fluctuations case and exponential correlated case 



with the smallest correlation radius R c — 5A, for both 
initial states and temperatures. And the dimerized cor- 
relations result in efficiencies that are similar to those 
produced by the exponential correlation with R c = 10A. 
This is not surprising since the dimers in FMO are formed 
by pigments separated by roughly 10A. Also, in compar- 
ing Figs, [I] and [3] we see that energy transfer from initial 
state 6 is more sensitive to spatiotemporal correlations 
than from initial state 1. 

Comparison of Fig. [TJwith Fig. [2] and Fig. [3] with Fig. 
[4] shows that the effect of increasing temperature, which 
in our model increases the variance of fluctuations, is to 
render the excitation transport less sensitive to temporal 
and spatial correlations. The optimal temporal correla- 
tion time is less pronounced; there is a wide plateau of 
comparable efficiencies across the range 30fs < r c < 90fs 
(Figs. [2] and [4]). Similarly, the variation of average trap- 
ping probability with spatial correlation is less at high 
temperature, although the general trend of decreased ef- 
ficiency with increased spatial correlation still persists. 
These simulations also show that the average efficiency 
of transport is more robust to the initial state at higher 
temperatures - i.e. there is less variation between Figs. 
[2] and [4] than between [T] and [3] 

For completeness, in Fig. [5]we also show representative 
time traces of the average trapping probability, for tem- 
poral correlation times r c = 45 fs and t c — 120fs, and 
for the three different spatial correlation models. The 
trapping probability curves have similar trends for all 
spatiotemporal correlations and the insets show that the 
dependence on spatial correlation length becomes less 
pronounced at higher temperatures. Fig. [6] shows the 
survival probability: P surv = 1 — Ptrap at the two tem- 
peratures 77K and 300K, for the case of no spatial cor- 
relations and temporal correlation time r c = 45fs. The 
time dependence of this survival probability is often ex- 
amined in the analysis of random walk dynamics with 
trapping [33]. For the stochastic Liouville description 
used here and for the graph defined by the FMO com- 
plex, we see that at T = 77K the survival probability at 
short times shows a complex modified exponential decay 
that is well fit (except in the t — ¥ limit) by a function 
of the form: P surv oc e^ ;Ci * with i ranging from to 3 
at least, containing both subexponential (cj < 0) and su- 
perexponential (c, > 0) contributions. The decay at long 
times fits well to a simple exponential decay e - ^ 4 . At the 
higher temperature T = 300K, the modified exponential 
decay regime is shortened significantly, and most of the 
decay is simply exponential. 



III. DISCUSSION 

To understand the behavior of trapping efficiency with 
temporal correlation, we focus on the case of no spatial 
correlation in fluctuations since the general behavior with 
respect to r c is the same for all spatial correlations con- 
sidered. We begin by examining the excitonic structure 
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FIG. 3: (Color online) Trapping probability at 20 ps as a function of the temporal correlation of the FMO 
environment, for various spatial correlation models. The initial state is an excitation localized on BChla 6, and the 
temperature in both panels is T=77K. The curves are polynomial fits to the data points that indicate the general 
trend, and the error bars show the standard deviation of the average taken over 100 sample evolutions. The color 
coding for the curves is the same as in Fig. [TJ Panel (a) compares the trapping probability for the case of no spatial 
correlation to the case of dimerized correlations, and panel (b) compares no correlation with exponential correlations. 
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FIG. 4: (Color online) Trapping probability at 20 ps as a function of the temporal correlation of the FMO 
environment, for various spatial correlation models. The initial state is an excitation localized on BChla 6, and the 
temperature in both panels is T=300K. The curves are polynomial fits to the data points that indicate the general 
trend, and the error bars show the standard deviation of the average taken over 100 sample evolutions. The color 
coding for the curves is the same as in Fig. [TJ Panel (a) compares the trapping probability for the case of no spatial 
correlation to the case of dimerized correlations, and panel (b) compares no correlation with exponential correlations. 



of FMO, shown in Table [TJ The wavefunctions of the two 
lowest energy excitons have significant overlap with the 
trapping site, BChla 3. Similarly, the initial state of an 
excitation localized on site 1 (site 6) primarily has com- 
ponents of excitons 2 and 5 (excitons 4 and 6). Using 
a second-order cumulant expansion technique [48], Blu- 
men and Silbey have shown that for such a model, when 
written in the interaction picture with respect to H e \, the 
populations in the exciton basis evolve according to the 



master equation 



dt 



Paa(t) 



1 



N-l 
,8=0 



where N is the total number of chromophores, and we 
have used Greek indices for the exciton basis. r Q( g is a 
population transition rate from exciton level a to ft, and 
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FIG. 5: (Color online) Average trapping probability as a function of time for calculations with two representative 
environmental correlation times, (a) t c = 45fs and (b) r c = 120fs. In each case the main panel shows the behavior at 
T = 77K while the insets show behavior at T — 300K. The initial state is an excitation on BChla 6. In both panels 

the blue (solid) curves are for no spatial correlations, C", and the red (dotted) and black (dashed) curves are for 
exponential spatial correlations with R c = 10A and R c = 20A, respectively. The averages are taken over 100 sample 
evolutions. Error bars are omitted for clarity: the variation from these average curves is small. 




is explicitly (in our notation): 



FIG. 6: (Color online) Log-linear plot of the average 
survival probability as a function of time for 
calculations with environmental correlation time 
t c = 45fs and initial state an excitation on BChla 6. 
This figure only shows the case of no spatial 
correlations; the temporal scaling is similar for all three 

models of spatial correlations analyzed in this work. 
The blue (solid) curve shows behavior at T = 77K and 
the red (dotted) for T = 300K . The inset is a zoom to 
the short time region. The average is taken over 100 
sample evolutions. Error bars are omitted for clarity: 
the variation from these average curves is small. 



2A 



1 + Tc 2 <8 



uj aj3 sin(w a/3 t) 



cos(a> Q ^i) 



where uj a p = uj a — up is a difference between exciton 
frequencies. For times t 3> r c , which is when most of 
the excitation transfer occurs, the second, time depen- 
dent term in Eq. ([2| is damped by the exponential pref- 
actor and can be neglected. Hence the transition rates 
are essentially determined by the constant term, which 
is equal to the Lorentzian power spectrum of the expo- 
nentially correlated noise: J\{ui) = ppp^ ■ In order 
to maximize transport efficiency it is advantageous to 
maximize the rate of transitions between exciton states. 
The time-independent part of the rates T a p is maximized 
when the correlation time of the fluctuations matches the 
energy differences between the exciton levels of the com- 
plex - i.e., t c = l/ujafj. Physically, this results from 
the fact that noise power at an exciton energy difference 
drives population transitions between the corresponding 
exciton levels. As can be seen from Table ffl the exci- 
ton energy differences in FMO are mostly in the region 
AE ~ 90 — 350cm -1 . For this range of energy differ- 
ences, the range of correlation times that maximizes Ja 
is t c ~ 15 — 60fs. This is precisely the range in which 
the peaks in average efficiency lie in Figs. [l]-|4j Thus we 
conclude that the peak in average efficiency with respect 
to temporal correlation is due to noise assisted transfer. 
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(a) T = 77K (b) T = 300K 



FIG. 7: (Color online) Average total coherence (main axes) and average displacement in A (insets) in the FMO 
complex as a function of time for two temperatures. The initial state is an excitation on BChla 6 and r c = 45fs. The 
blue (solid) curves are for no spatial correlations, C^; the red (dotted) and black (dashed) curves are for 
exponential spatial correlations with R c = 10A and R c — 20A, respectively. The average is taken over 100 sample 
evolutions. Error bars are omitted for clarity, the variation from these average curves is small. The inset shows the 
mean displacement X (defined in the main text) as a function of time for the same spatial correlation cases. 
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2 
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Energy 





102.8 


177 


272.7 


297.5 


402.7 


497.2 


Site overlap 


3,4 


3, 4, 5, 7 


1, 2 


5, 7 


4, 5, 6 


1, 2 


5, 6 



TABLE I: FMO excitons (energy eigenstates of H e i). 
The exciton energies are in cm" 1 (normalized so that 
the lowest energy exciton has energy zero). The third 
row shows the sites that the exciton wavefunction has 
largest overlap with. 



Now we turn to the effects of spatial correlations. The 
primary mechanism by which energy fluctuations affect 
excitation transport is by modulating energy mismatches 
between chromophores [ITJ [19]. The rate of excitation 
transfer between two chromophores is enhanced when 
their average energy gap decreases. Positive correlations 
in the energy fluctuations suppress line broadening and 
in the limit of perfectly correlated fluctuations the en- 
ergy gaps remain unchanged. This can equivalently be 
viewed as a renormalization in which positive correlations 
reduce the reorganization energy, and hence the dephas- 
ing rate |34j . The suppression of line broadening has two 
key effects on excitations: it leads to longer-lived coher- 
ence between sites and it slows the average transport by 
reducing transfer rates between energy mismatched chro- 
mophores. We illustrate both of these effects in Fig. [7] 
The main axes plot the time evolution of total coherence 
in FMO, which is defined here as: C — \Pij\- The 

initial state is an excitation localized on BChla 6, and 
the two panels show the decay of coherence at two tem- 
peratures. Both plots clearly show the preservation of 
coherence by spatially correlated fluctuations, although 



the effect is less dramatic for the higher temperature. 
Further, to demonstrate the second point that correlated 
fluctuations lead to slower excitation transport, in the 
insets we plot the mean displacement of the excitation, 
defined as X{t) = J2i^6iPii{t), where dei is the distance 
between the initial site 6 and site i. 

It is clear from these plots that after the first few hun- 
dred femtoseconds, the average rate of transfer of the 
excitation is reduced in the presence of spatially corre- 
lated fluctuations. Again, the effect is less dramatic at 
higher temperatures where the line broadening is inher- 
ently larger and the renormalization by spatial correla- 
tions has less of an impact. Thus we conclude that pos- 
itive spatial correlations in environmental noise lead to 
both longer-lived coherence and slower excitation trans- 
port. These conclusions agree with recent studies of the 
influence of spatial correlations in chromophoric systems 
using a variety of techniques [9[ EH ES EEl [39] . 

As noted above, in view of the fact that the exact na- 
ture of the spatial correlation present in the environ- 
mental fluctuations is unknown [THl HD], we also sim- 
ulated excitation dynamics assuming a polynomial in- 
stead of exponential dependence of the spatial correla- 



tions. In particular, we used C 



DB 



l 

d 2 . 



X 0, 



i.e., an 



inverse squared dependency of vibrational correlations 
in the protein scaffolding on the distance between chro- 
mophores. Here f3 is the largest constant < 1 that ensures 
that the matrix C DB is positive definite: the magnitude 
of j3 thus constitutes an upper bound on the magnitude of 
correlations between fluctuations on different sites. With 
the inter-chromophore distances from Ref. [23], we find 
(3 ss 0.85. Our calculations of the dynamics resulting 
from this correlation of the environmental fluctuations 
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(not shown here) are very similar to that presented for 
exponentially decaying correlations in Figures 1-5. In 
particular, there is an optimal temporal correlation time 
of ~ 40fs for initial state 1 (and ~ 30fs for initial state 
6), and the overall efficiency of energy transfer is seen to 
be lowered by the spatial correlation, while coherence in 
the site basis is enhanced, as in Figure [7] 

IV. CONCLUSION 

We have systematically studied the effects of temporal 
and spatial correlations of noise on the transport of ex- 
citations across a prototypical light harvesting complex. 
We emphasize that since we employed a classical stochas- 
tic model of the phonon environment, it is not possible 
to draw quantitative conclusions regarding energy trans- 
fer in FMO from the present study. However, our model 
does enable us to examine the relative effects on trans- 
port of varying amounts of correlation. We have shown 
that temporal correlations can enhance noise power in 
certain spectral regions and consequently enhance exci- 
ton transitions lying in these regions. We also showed 
that spatially correlated fluctuations can preserve coher- 
ence while at the same time resulting in slower transport 
to the trapping site. While this paper was under review 
we learned of recent studies using the more physically 
accurate generalized Bloch-Rcdfield equation approach 
[351 El] that arrive at similar conclusions. These authors 
further investigated the dependence of these phenomena 
on the reorganization energy of the protein pigment com- 
plex and found that depending on the magnitude of the 
reorganization energy, positive spatial correlations can 
either decrease or increase the overall efficiency of the 
energy transfer process. This interesting result is con- 



sistent with both the fact that spatial correlations act 
to effectively decrease the reorganization energy of the 
system and the fact that there is an optimal value for 
the reorganization energy. Thus depending whether the 
reorganization energy is smaller or larger than the op- 
timal value, spatial correlations may shift the value ei- 
ther towards the optimal, increasing the efficiency, or 
away from it, decreasing the efficiency. Finally, although 
we have specifically analyzed the FMO complex in this 
work, the conclusions drawn here about the effects of 
correlated fluctuations on excitation transport will also 
apply to larger light harvesting complexes and are gener- 
ally applicable to transport phenomena in densely packed 
molecular assemblies, including J-aggregates and other 
photosynthetic units [37] . 

This study also raises several intriguing questions of 
biological import. Most significantly, why is it that spa- 
tially correlated fluctuations are likely present in some 
LHCs [1, 2 J if they serve to reduce rather than to enhance 
the efficiency and speed of excitation transport? Second, 
is the correlation time of environmental fluctuations in 
FMO within the window that maximizes transport effi- 
ciency? Clearly, more detailed studies are warranted of 
the complex environments for electronic energy transport 
that are found in natural organic molecular assemblies, 
in particular for light harvesting complexes. 
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